Structural behavior of uranium dioxide under pressure by LSDA+U calculations 
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Structural behavior of UO2 under high-pressure up to 300 GPa has been studied by first-principles 
calculations with LSDA+U approximation. The results show that a pressure induced structural 
5^ ■ transition to the cotunnite-type (orthorhombic Pnma) phase occurs at 38 GPa. It agrees well with 

experiment observed ~42 GPa. A new iso-structural transition following that is also predicted taking 
place from 80 to 130 GPa, which has not yet been observed in experiments. Further high compression 
beyond 226 GPa will result in a metallic and paramagnetic transition. It corresponds to a volume 
' of 90 A 3 per cell, in a good agreement with previous theoretical analysis in the reduction of volume 

required to delocalize 5/ states. 
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Uranium dioxide (UO2) is widely used as fuel material in nuclear reactors. It was extensively studied from 1940th 
in physical, chemical and thermodynamic properties by experimental and theoretical methods (especially during last 
decade) due to its important applications^i 3 -^^ The shortage of energy around the world makes us ask for more 
contribution from nuclear power, where the burn-up efficiency of nuclear fuels is a bottle-neck. In order to tackle 
this difficulty, understanding the detailed behavior of fuel materials under burning-up and irradiation is important. 
Q ■ Recent use of fuel materials to high burn-up shows many microstructure formations, which is not possible to access 
O \ by empirical approach and atomic scale theoretical analysis is highly requested!^ 

Previous theoretical studies on UO2 mainly focused on defects effects arised from irradiation damages 6 ' 9 and ther- 
modynamic properties near ambient pressur e 1 ' 10 ! 11 with semi-empirical approaches. Electronic properties and large 
, scale intrinsic structural behavior under disturbance (say, compression, tension and distortion of lattice) were rarely in- 
' vestigated in despite of its importance in most properties of UO2. This partly is due to the lacking of a reliable method 
\ to deal with this kind of complex materials before. The development of density functional theory (DFT) changed the 
situation greatly and provides a quantum mechanics based theoretical approach to tackle this problem! 1 1 

However, to our knowledge only a few ab initio electronic structure studies have been published on UO2, most 
of which were based on conventional LDA or GGA approximation of the exchange-correlation energy i 3 ' 4 ' 15 i 16 i 17 i 18 It 
is well recognized that strong Coulomb correlation among partly filled / electrons of uranium atoms makes these 
approximations failed. Usually a metallic ground state is predicted for UO2 instead of the experimental observed 
antiferromagnetic semiconductor i 19 i 20 The same problem exists for transition metal oxides, and raises questions about 
the applicability of DFT approach to these materials. Fortunately, a method combining spin-polarized local density 
approximation (LSDA) and on-site Coulomb repulsion among localized d or / electrons / 3 ' 14 namely, LSDA+U method 
was proposed and has shown its capability to treat this problem ! 21 ' 22 

Usually a Hubbard Hamiltonian with two empirical parameters is employed to describe the Coulomb interaction 
between 5/ electrons localized on uranium sites in UO2. Adding this Hamiltonian to the conventional LSDA (or CCS) 
energy functional, one arrives at a point where all orbitals except those included in the model Hamiltonian are treated 
within the framework of LSDA (or GGS) while the localized 5/ states are treated by the unrestricted Hartree-Fock 
(UHF) approximation, 22 namely, 
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E L sda + u = E LSDA [{ Si }\ + V^£L p%p% (1) 

1,3," 

where is the density matrix of electrons occupying a partly filled electron shell (5/ in UO2), a refers to spin 
direction, and {ei} is the Kohn-Sham eigenvalues. Self-consistent solution of DFT with this energy functional gives 
that strong correlation effects associated with 5/ states are going to affect all other states as well, in particular though 
2p states of oxygen do not influenced by Hubbard correlations directly, they are really linked to localized 5/ states via 
hybridization terms. It is necessary to point out that LSDA+U method is not a self-determined approach. The results 
depend on model parameters U and J very much, which should be chosen carefully by comparing with experimental 
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data. Fortunately this can be done very well with just a small set of data and preserve the predicability of the method 
mostly^iS 1 .^ 

There have been several works with LSDA+U approximation on uranium dioxide published. All of them were 
near the equilibrium volume at ambiance pressure for fluorite structure and focused mainly on electron energy loss 
spectr a 21 ! 22 ' 23 and magnetic structured These calculations showed that the results of LSDA+U in a well agreement 
with experiments. However, no attempt was made to investigate the structural behavior of UO2 under pressures with 
LSDA+U method, which may be fundamental for understanding the behavior of nuclear fuel under irradiations. By 
far the validity of LSDA+U method beyond fluorite structure for UO2 has not been confirmed yet. Recent hydrostatic 
compression experiment 2 ^, makes it possible to check it by comparing with measured equation of state. On the other 
hand, first principles calculations without Hubbard correction on GGA(S) approximation showed that it can give 
almost correct energy information for UO^^ regardless a wrong electronic band structure was predicted. Specially, 
by calculating the lattice parameter and bulk modulus of fluorite structure UO2 with various approximations, J. 
C. Boettger argued that density gradient corrections, spin polarization and spin-orbit coupling effects are equally 
important, and suggested when only structural properties are concerned LSDA+U is not necessary.— However the 
predicted wrong ferromagnetic ground state weakened the creditability of his argument. Other calculations ignored 
spin-orbit coupling also gave reasonable lattice parameter and bulk modulus) 3 -^ 8 - indicates spin-orbit coupling is not 
so important for this case (though a large impact on magnetic property is expected). We will show in this paper 
that it should be careful when GGA(S) approximation is used because the coincidence of cohesive energies of UO2 
calculated by GGS with LSDA+U approximations is valid only for fluorite phase. An energy difference will appear if 
other structures are involved. 

In this paper we will study the structural stability of fluorite phase (with Fm3m space group) and cotunnite phase 
(Prima space group) of uranium dioxide under hydrostatic pressures using DFT method based on LSDA and GGS 
approximations plus Hubbard correction. Calculation methodology is presented briefly in next section. We will discuss 
a little bit about the widely used rule of common tangent of energy curves to determine the transition pressure of 
pressure-induced structural transition, because this rule fails in a case when an energy barrier existed. A more general 
rule is proposed, which can give the energy barrier when experimental transition pressure is available. Finally, a 
detailed comparison of our results with static high-pressure experiments is given, associating with a discussion on 
ultra-high pressure behavior of UO2 crystal. 



II. METHODOLOGY 



Total energy curves of both phases (Prima and Fm3m) at different volumes are computed with VASP code i 26 ' 27 The 
Prima structure is fully relaxed to get all Hellman-Feynman forces smaller than 0.002 eV/A, while fluorite structure 
keeps the ideal geometry due to all coordinates are completely determined by the symmetry. For comparison, both 
spin-polarized generalized gradient approximation (GGS) 28 and local density approximation (LSDA) 2 ^ with/without 
Hubbard U term energy functional are used. The parameters of Hubbard term are taken as U = 4.5 eV and J = 
0.51 eV, which was checked carefully by S. L. Dudarev et al. for fluorite U02j 2 ° i21 ' 22 Calculations employ projector- 
augmented wave (PAW) pseudopotential a 30 ' 31 with a cutoff kinetic energy for planewaves of 400 eV. Integrations in 
reciprocal space are performed in the first Brillouin zone with 18 irreducible k-points for fluorite structure and at 
least 28 irreducible k-points for cotunnite phase generated with the Monkhorst-Pacb 3 ^ scheme. Its convergence is well 
checked. The energy tolerance for the charge self-consistency convergence is set to 1 x 10~ 5 eV for all calculations. 
Cohesive energies at different volumes are extracted from the total energies by subtracting spin-polarized isolated 
atom contributions. Then, they are fitted to a Morse-type energy function 



E(V) = D 



-|[(v/i/ ) 1/3 -i 



1 -1 



(2) 



to facilitate post-analysis. It is necessary to point out that for Prima phase we also used a different U value obtained 
by fitting to experimental data of Prima phase since U — 4.5 eV fails to predict the correct transition pressure. This 
implies that structure or lattice distortions would have considerable impact on on-site coulomb interaction. For the 
same structure, however, we find the dependence of U on pressure is ignorable. 

The equation of state (or compression curve) at zero-Kelvin is calculated directly by an infinitesimal variation of 
cohesive energy with respect to volume given by P — —dE/dV. Usually the phase transition pressure is determined 
by the common tangent of their energy curves, which can be derived simply as follows. At thermodynamic equilibrium 
state under finite pressure, the enthalpy must be minimized, i.e., SH = 0. In a case two phases in equilibrium, there 
is a variation of enthalpy with respect to the concentration of each phase besides with respect to volumes. The latter 
gives Pi = —dEi/dV (where i is phase label) and the former results in P = —AE/AV — (E2 — Ei)/(Vi — V2) with 
5H = 8x(AE + PAV), where Sx is the concentration variation of, say, the first phase and AE (AV) is the energy 
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TABLE I: Cohesive energies of uranium dioxide at GPa. 
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(volume) difference between these two phases. The balance condition of pressure requires P — Pi (for i — 1, 2) 
namely, 



dEj 
dV 



v=v t Vi - V, 



E ^ El (3) 



It is exactly the common tangent rule for transition pressure of pressure induced structural transitions. Evidently, 
PAV provides the least energy AE required to drive a transition from phase 1 to phase 2. The transition pressure 
equals to P if no energy barrier exists, which is a common case for usual metals and alloys. However, when an energy 
barrier with an amplitude of Aw is involved ) 33 ! 34 ' 35 the work done by external pressure P' should be large enough to 
get over the barrier in addition to the energy difference AE. Then the variation of enthalpy with respect to phase 
concentration should be SH = (AE + Aw + P'AV) ■ Sx = 0. Obviously the common tangent rule becomes invalid 
here. The hysteresis pressure is given by AP = P' — P. Without knowledge about the energy barrier, one cannot 
determine the transition pressure P' by energy curves itself. However, in contrast, one can deduce the energy barrier 
amplitude with measured transition pressure P' by 

Aw = -(AE + P'AV). (4) 



III. CALCULATIONS AND DISCUSSIONS 



A. Cohesive energy 

Calculated cohesive energies with different approximations for fluorite (Fmim) and cotunnite (Pnma) structures 
of uranium dioxide as well as the parameters fitting to Eq.([2]) are listed in table [H where the energy is for per atom 
and equilibrium cell volume is given by Vq = r 3 . The cohesive energy for a cell of U4O8 is given by multiplying D with 
12. For comparison purpose, other calculated and observed value a 17 i 18 ' 21 i 25 ' 36 are also listed. It should be noticed 
that Pnma phase has a smaller effective cubic lattice constant and bulk modulus than fluorite phase at zero-pressure, 
which implies it will become stable under compression. 

Variation of cohesive energy of Fm3m phase along cubic lattice constant (V 1 ^ 3 ) is shown in figure[TJ It is interesting 
to see that GGS and LSDA+U give quite similar energy curves, confirming previous calculations that GGS also can give 
reasonable energy information for fluorite phase of UO2 in spite of the corresponding electronic density of state (DOS) 
is wrong4 i 15 ' 18 i 37 However, we should emphasize here that it is just a coincidence. Analogous to the case without 
Hubbard correction, LSDA+U overestimates the binding energy slightly and GGS+U underestimates it. On the other 
hand U term uplifts the binding energy wholly, and results in this coincidence. We can also see from table fl] that 
PAW method outperforms ordinary pseudo-potentials both for GGS and LSDA approximations in terms of equilibrium 
volume, cohesive energy and bulk modulus ^ 17 ' 18 Furthermore, Our LSDA+U calculations with PAW potentials give 
results in perfect agreement with experiment s 21 ' 25 ! 36 (in particular the calculated equilibrium lattice constant of 5.44 A 
vs observed 5.46 A and bulk modulus of 208.3 GPa vs 208.9 GPa). It also predicts an antiferromagnetic ground state 
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FIG. 1: Comparison of cohesive energy curves for fluorite structure of UO2 calculated by LSDA, GGS, LSDA+U and GGS+U 
approximations, respectively. Notice GGS and LSDA+U approximations give very close energy, especially at high-compression 
region. 




Cell volume (A 3 ) 

FIG. 2: Comparison of cohesive energies of Pnma and Fm3m phases along cell volume. A phase transition at 7.8 GPa is 
predicted by the slope of common tangent rule for (7=4.5 eV CeLSG, ELS shown in the inset. 

with a band gap of ~1.45eV, agrees with previous calculation very well^ To reproduce the X-ray photoelectron 
spectroscopy^ observed band gap of ~2 eV, Dudarev et al. argued that to take spin-orbit coupling into account^ 
is necessary. We confirmed this by a spin-orbit coupling calculation implemented in VASP which gives a band gap 
of 2.04 eV. GGS+U approximation, however, gives a larger equilibrium lattice constant and smaller bulk modulus, 
despite the cohesive energy is more close to the observed value, as well as a band gap of 1.6 eV. Totally speaking, 
LSDA+U outperforms GGS+U approximation for this set of U term parameters. It is necessary to point out that the 
discrepancy with previous LMTO calculation 21 should be owing to their convergence precision is not so good. Their 
calculation gave quite poor mechanical properties^ that implying the force is inaccurate. Later calculation by the 
same authors improved this— The spin-orbit coupling is ignored in our following calculations. The resulting error 
can be estimated at a lattice constant of 5.44 A for fluorite phase, where spin-orbit coupling decreases the cohesive 
energy about 0.3 eV for per atom, very close to the energy difference between GGS and LSDA+U at the same volume. 
Thus we can expect that neglect spin-orbit coupling will gives an error of 0.012 A in lattice constant and 3 GPa in 
bulk modulus, the same difference as GGS. 
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FIG. 3: Calculated compression curves of uranium dioxide along relative volume compared with experimental measurements. 
An unexpected large discrepancy of transition pressure between measurement (point B) and theoretical prediction (point A) 
is obtained for U45 case. However, a better result is recovered with £7=6.0 eV. 



B. Cotunnite phase 



To optimize the geometry of Prima phase at different pressure, LSDA+U method with U = 4.5 eV is employed. To 
avoid the Pulay stress problem (which arised from the fact that the plane wave basis set is not complete with respect 
to changes of the volume), structure relaxation calculations are performed at fixed volumes rather than under constant 
pressures. Then pressure is derived from the energy-volume relation. The structure is fully relaxed to optimize all 
internal coordinates and cell shape, while the symmetry of Prima space group is kept. Calculated cohesive energy 
curve is shown in figured] For comparison the curve of fluorite phase is also given as dash-dotted line. It shows that 
under high pressure Prima phase becomes stable. A transition pressure of 7.8 GPa is given by the slope of common 
tangent as showing in the inset. This value is quite unexpected because it is less than 1/5 of the experiment observation 
as ~42 GPa.— It is very small even if compared with another early measurement that predicted a pressure-induced 
phase transition to orthorhombic Cmcm phase at ~29GPa— i (which has not yet been repeated by other authors). 
Nevertheless, the calculated volume reduction of 6.4% agrees well with the observed 7% at the beginning of cotunnite 
phase.— 

Then one may ask what is the matter with it? Is the LSDA+U approximation wrong? From table U and the 
comparison of its results at equilibrium volume with experimental data for fluorite phase, we do not think so. Actually, 
as figure [3] shows, LSDA+U gives a compression curve that agrees very well with experiment o 25 ' 41 for Fmim phase, 
which means that the Hubbard £7 parameter is reasonable and insensitive to pressure. Clearly we cannot attribute this 
deviation to the failure of density functional theory or LSDA+U approximation. Figure [3] also shows the P-V curve 
calculated with GGS approximation. It is worse than LSDA+U and the transition pressure is also as low as 32 GPa. A 
hysteresis pressure about 34 GPa is estimated by using the transition pressures observed in experiment and calculated 
with U = 4.5 eV, which are marked by arrows B and A in figure [31 respectively. As discussed in previous section, this 
hysteresis of transition pressure would imply an energy barrier existing. In fact it is very common for ionic crystal 
and semiconductors. 33 ' 34 ' 35 For example a phase transition of GaN from wurtzite to rocksalt phase, where a large 
hysteresis of pressure is observed. By using Eq.(|4]), the cohesive energy curves of Prima and Fm3m phases (Eq.([2]) 
and parameters listed in table [I]), and the experimental transition pressure of 42 GPa? 2 ^ we estimate an energy barrier 
as ~2.1 eV per cell (U^Os)- This value is large enough to survive Prima phase to ambient condition. Unfortunately, 
no experiment shows this event. 

Therefore, the only possibility is that the discrepancy in Prima phase resulted from the dependence of £7 on structure 
(or lattice distortions). We obtained a different £7=6.0eV by fitting to the measured P-V data of Prima phase. The 
resulting energy curve and Morse function parameters are given in figure [21 and table HI respectively. We can see 
U = 6.0 eV gives a quite similar energy curve as U = 4.5 eV, except the wholly uplifting of the curve. Hereafter, 
all calculations will be performed for £7=6.0 eV and £7=4.5 eV separately, so we assign the former case as £76 and 
the latter as £7 45 for briefness. Although the improvement on P-V curve in £76 is limited, as figure [3] shows, the 
calculated transition pressure is corrected to ~38 GPa, almost five times of U 45 case and in a good agreement with 
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FIG. 4: Electronic density of state for 5/ states of UO2 at a cell volume of 131.4 A 3 . Transition to Pnma will increase the 
band gap and shrink the energy region of localized states. 




FIG. 5: Variation of magnetic moment of uranium atoms with cell volume for Pnma and Fm3m phases. The metallic transition 
is indicated by the arrows, where a paramagnetic transition also occurs simultaneously. 



observed 42GPa. The resulted hysteresis pressure is just 4GPa, which ends up an energy barrier as 0.018 eV/atom 
and ignorable at room temperature. Obviously U6 is more credible than U 45 since it is compatible with the fact 
that no Pnma phase has been observed under ambient condition. The calculated reduction of volume at transition 
from Fm3m to Pnma phase is 6.2%, close to {745 case, also agrees well with experimental data. 

Figure Q] compares the density of state (DOS) of 5/ states in Fm3m and Pnma phases of UO2 at a cell volume of 
131.4 A 3 , close to the transition pressure of U6 case. The most remarkable difference is the increase of band gap from 
0.8 eV in fmorite phase to 2.4 eV in Pnma phase. As a consequence, unoccupied states also move outwards. Below the 
Fermi level, different from fluorite phase where a nearly dispersionless band containing two well-localized 5/ electrons 
that lies roughly from -1.8 to OeV, in Pnma phase these localized states are further narrowed to start from -1.0 eV, 
while the valence 5/ state is expanded from -3.7 to -1.0 eV, too. To completely delocalize the localized 5/ states, a 
pressure above 121 GPa is required for C/45 and beyond 226 GPa for U6, we will discuss this in next subsection. 
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FIG. 6: Total electronic density of state calculated with 17=6.0 eV for Prima phase under a pressure of 226 and 294 GPa, 
respectively. Here the Fermi level is at eV. 



C. High pressure behavior 



The variation of local magnetic moment of uranium atoms with cell volume is almost the same for Pnma and 
FmSm phases in Z745 case, implying the magnetic property is insensitive to structural transition in UO2. As figure [5] 
shows, in despite of that GGS+U and LSDA+U approximations give much different cohesive energies for the fiuorite 
phase, the calculated magnetic moment of uranium atoms is very close for a large range of volume, except for the 
highly expansion region (V > 300 A 3 ) where atoms trend to be isolated. At equilibrium volume, our calculation gives 
a moment of ~1.93^s in a good agreement with previous calculation^! and slightly larger than observed 1.74 ^3, 42 
It can be improved by including spin-orbit coupling to 1.88 ub with an orbit contribution of 0.46 ub- This value is 
much smaller than all-electron calculation where an orbit moment of 3.6 ub was predicted^ 

As shown in figure[5l there is a flat level for the local magnetic moment of uranium within moderate pressure range. 
A transition from antiferromagnetism to paramagnetism is observed at a volume between 102. 1^108. 4A 3 for t/45 
case (equivalent to 121 and 159 GPa in pressure). It corresponds to a volume of 63~66.4% of the equilibrium volume 
of fiuorite phase and an effective cubic lattice constant as 86^88% of the latter phase. Increase U to 6.0 eV postpones 
the paramagnetic transition to higher pressure as 226~294GPa, which has an effective cubic lattice constant that is 
82^84% of the fiuorite phase at ambient condition. It is worthwhile to point out that at the same volume a metallic 
transition also occurs due to completely derealization of 5/ states. Figure [H] shows the total DOS of Pnma UO2 under 
high pressures. We can see the band gap disappear completely between 226~294GPa. The transition volume is in a 
good agreement with previous intuitive analysis that a reduction in the effective cubic lattice parameter to 82% of the 
equilibrium lattice parameter do (of Fm3m phase) is required to have 5/ states in the conduction band 4 s - It is clear 
that the paramagnetic transition is driven by the derealization of the two pre-localized 5/ electrons which become 
itinerant at this volume, and it is quite reasonable that the metallization is always accompanied by a paramagnetic 
transition for materials analogous to UO2 where both band gap and local magnetic moment are attributed to the 
same localized states. 

Below the metallic transition, we also find a new iso-structural transition occurring between 80~130GPa for Pnma 
phase. Figure [7] shows the variation of relative lattice parameters of Pnma phase starting from respective transition 
pressure of t/45 and U6 cases. Drastic variations were observed for all lattice parameters between 75^121 GPa for 
U45 and 80~130GPa for f/6, where the smallest axis b has a strong rebound and the middle a is collapsed. At 
higher pressure, the variations of relative lattice parameters become smooth and approaches isotropic compression. 
It is a typical structural transition. For U 45 case one may wonder whether there is some relevancy between this 
transition and the metallic one because they adjoin closely in pressure. However calculation with U 6 shows that 
they are irrelevant. By the way, At low pressure the calculated variation of relative lattice parameters is different 
from experimental observation, where the smallest axis b is most compressible whereas the a axis is most rigid. We 
do not know the exact reason for this discrepancy at present. But the experiment observed trend of relative lattice 
parameters cannot hold to high pressures because a stronger repulsive force will present along the shorter axis due to 
higher compression of electronic states. One can expect a rebound of the smallest axis at higher pressure. 
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Pressure (GPa) 

FIG. 7: Pressure behavior of relative lattice parameters of Pnma phase, where the drastic change in relative lattice constants 
(region between dotted lines) indicates an iso-structural transition. The curves with open symbols are calculated with U =4.5 eV, 
and those with filled symbols are for (7=6.0 eV. 




Intermediate structure index 

FIG. 8: Relative energies of intermediate structures interpolating Fm3m and Pnma phase linearly. The solid line without 
symbol is obtained with a linear interpolated value of U between 4.5 and 6.0 eV successively. 

D. Intermediate structures 

As discussed in previous subsections, The value of U depends on structure. This raises a question about the 
applicability of LSDA+U method to intermediate process of structural transition, since the energy is affected by this 
term directly. It is impossible to fit the U value for all intermediate structures with experimental data. Therefore, 
if one attempts to approximately model the transition (or lattice distortions) with just single or several values of U, 
evaluating the corresponding error becomes important. We do this job for UO2 by calculating the energy variation 
along the linear interpolated intermediate structures between Fm3m and Pnma phases under ~8 GPa, namely, a 
candidate transition path for U 45 case. In this calculation, no structure optimizing is performed. 

The result is shown in figure[8l where the respective energy of Fm3m phase is set as reference point for GGS, LSDA 
and a classical pair potential model.— For LSDA+U, only the energy of fluorite phase calculated with [7=4.5 eV is 
set as reference energy, to take varying U effect into account. Since U6 fails to model Fm3m phase and f/45 fails to 
describe Pnma phase, as the first level approximation, we interpolate the value of U between these two phases linearly. 
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The result is given in figure [8] as the solid curve without symbol. As expected, U 6 performs well for intermediate 
structures near Prima phase while f/45 becomes better for those close to fluorite phase. The largest error is 0.15 and 
0.14eV per atom for U6 and f/45, respectively. What amazing is that the classical pair potential model outperforms 
GGS/LSDA approximation in this test. The former has an error as 0.18 eV per atom and the latter two are 0.32 and 
0.34 eV per atom, respectively. It is about two times larger than LSDA+U approximation. This result of GGS/LSDA 
approximation is somewhat disappointed. For an unit cell of U4O8, it would lead to an error about 4eV in cohesive 
or formation energy. In this sense, the point defect formation energy calculated by M. Freyss et al*& is inaccurate 
and need further improvement with LSDA+U method in uranium defects case due to the large structure distortions. 
Finally, we would like to point out that the previous conclusion made by J. C. Boettger that "the strong correlation 
effects that are generally believed to produce the observed band gap do not have a significant impact on the binding 
properties of UO2" 15 should be treated carefully depending on the studied structures and required precision. 



IV. CONCLUSION 



The structural behavior of uranium dioxide under pressure up to 300 GPa was investigated by DFT method with 
GGS/LSDA approximations plus (or not) Hubbard U correction for strong correlated on-site Coulomb interactions. 
Comparison with experiment showed that LSDA+U gives the best description for UO2 in fluorite phase. However 
the calculated transition pressure to Pnma phase with the same U parameter was quite low, indicating the value of 
U depends on structure or lattice distortions sensitively. A better value of U for Pnma phase is obtained, which 
removes the factitious energy barrier predicted by t/=4.5eV. The error due to varying of U is estimated as just half 
of the error given by GGS/LSDA approximation, showing LSDA+U is more reliable. Higher pressure leads to an 
iso-structural transition followed by a metallic-paramagnetic transition, which takes place between 226^294 GPa with 
an effective cubic lattice parameter as 82~84% of the fluorite phase's at zero pressure, in a good agreement with 
previous theoretical analysis. 
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